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Transport as well as magnetic relaxation properties of the mixed state were studied on strongly underdoped 
Yi-xPrxBa2Cui07-i crystals. We observed two correlated phenomena - a coupling transition and a transition to 
quantum creep. The distribution of transport current below the coupling transition is highly nonuniform, which 
facilitates quantum creep. We speculate that in the mixed state below the coupling transition, where dissipation 
is nonohmic, the current distribution may be unstable with respect to self-channeling resulting in the formation of 
very thin current-carrying layers. 



I. INTRODUCTION 

The majority of experiments devoted to the study of 
dissipation in the mixed state of cuprate superconduc- 
tors were performed on optimally doped (maximum crit- 
ical temperature Tc) single crystals, thin films, or super- 
latttices. Meanwhile, in underdoped systems, the nor- 
mal state properties exhibit a number of unusual features 
which quite possibly hold a key to the understanding of 
the nature of the normal and superconducting states in 
cuprates. Also, as a result of the lower T^ and upper crit- 
ical field Hc2 , the dissipation in underdoped systems can 
be measured down to much lower reduced temperatures 
T/Tc in relatively small fields of a few tesla. In contrast, 
in optimally doped single crystals, the dissipation falls 
below typically detectable levels at much higher reduced 
temperatures. Thus, underdoped cuprates facilitate the 
study of a broader range of the magnetic field - temper- 
ature _ff — T phase diagram of the "vortex matter" than 
optimally doped superconductors. These facts provide a 
compelling reason to undertake a comprehensive study 
of the mixed state of strongly underdoped cuprates, in- 
cluding both transport and magnetic relaxation measure- 
ments. 

One of the outstanding and extensively debated ques- 
tions is the nature of the coupling transition in layered 
vortex systems. A sharp coupling transition has been 
observed in superconductor/insulator multilayers begin- 
ning with the pioneering work of Giaever M, and later 
in Refs. '^^- In a magnetic field H applied normally 
to the planes, the positions of 2D vortices (pancakes) on 
neighboring superconducting layers tend to be uncorre- 
lated at high temperatures. This vortex state is called 
2D liquid. As the temperature decreases, the correlation 
in the direction of the applied magnetic field strengthens 
and the vortices tend to form coherent 3D fiexible lines. 
For high- Tc superconductors, though, the existence of the 
coupling transition is not as obvious. Several groups us- 



ing a 6-terminal (flux transformer) technique have ar- 
rived at conflicting conclusions. 

Safar et al. |j] reported the observation of a transition 
to a 3D liquid in YBa2Cu^OT^s which manifests itself 
as a convergence of voltages generated by the motion of 
vortices on opposite faces of the sample. However, such 
a strong manifestation of the coupling transition appears 
to be the exception rather than the rule. Other groups 
have observed that voltages generated on opposite faces 
of Bi2Sr2CaCu20s |-|| and YBa2Cu^Oj^s [§ single 
crystals diverge with lowering temperature, rather than 
converge. 

A point of view which can reconcile the results of these 
experiments is that a signiflcant increase in the corre- 
lation length Lc of the vortices in the c-direction in- 
deed takes place at a well defined temperature, but Lc 
may remain smaller than the thickness of the sample. 
Thus, voltages generated on opposite faces of the sample 
may never converge in spite of a macroscopic correlation 
length. The underdoped cuprates present an opportu- 
nity to test this idea for reasons presented above. If the 
coupling transition can be proven in an underdoped sys- 
tem, it must also occur in optimally doped cuprates since 
the coupling between the Cu02 bilayers decreases with 
underdoping. 

A second important question in the physics of vortex 
matter is the possibility of non-activated, temperature 
independent, creep due to quantum rather than classical 
(over the barrier) relaxation at low temperatures. Super- 
conductors represent, perhaps, the only system in which 
relaxation due to quantum creep is an experimentally ac- 
cessible phenomenon. Here, a strongly non-equilibrium 
macroscopic metastable state relaxes coherently without 
thermal activation. In contrast, in the majority of other 
macroscopic metastable systems, relaxation proceeds as 
a sequence of a large number of uncorrelated microscopic 
steps, requiring thermal activation over an energy bar- 
rier. 

Yet, the evidence of a magnetic relaxation rate that 



does not extrapolate to zero as T ^ has not reached 
a point where experimental data can form a cohesive 
picture of the phenomenon. While non-vanishing mag- 
netic relaxation has been observed in both single crystals 
and thin films |lO-n3[ , non- vanishing resistance has been 
observed only in ultrathin films IIj-M. This has con- 
tributed to the assertion that temperature independent 
resistance in films and non-vanishing low-temperature 
magnetic relaxation in single crystals are unrelated phe- 
nomena. Since these two types of measurements involve 
different ranges of current (small currents in transport 
and large currents, close to the critical current Jc, in 
magnetic relaxation), they are open to alternative inter- 
pretations, not related to quantum creep. For example, 
Gerber and Franse ||lj] have argued that non-vanishing 
magnetic relaxation at low temperatures may result from 
self-heating, so that the local temperature of the sample 
is higher than that of the ambient. 

The best way to address these issues is to conduct both 
types of measurements on the same system. A signature 
of a temperature independent creep appearing in both 
transport and magnetic relaxation measurements at the 
same temperature would be convincing proof that this 
phenomenon is not an artifact and reflects a fundamental 
change in the relaxation process. 

The strongly underdoped Yi-a;PrxBa2Cu^OT-s (Tc ~ 
17—21 K) system is an excellent candidate for this study. 
One interesting aspect of this low Tc system is that its 
normal state, revealed by the suppression of supercon- 
ductivity by magnetic field, is a two-dimensional insu- 
lator pSl. The reduced dissipation in the normal core 
of a vortex (due to a large normal-state resistivity p„) 
increases the mobility of the vortices 11^ and, there- 
fore, favors quantum tunneling |2Q-p2]. Thus, in such 
an insulator-superconductor material we can expect that 
the transition from thermally activated to quantum creep 
takes place at higher temperatures than in conventional 
superconductors or in more metallic cuprates with higher 
Tc. 

In addition, as shown below, this insulator- 
superconductor system provides an example of a drastic 
departure from the current theoretical understanding of 
the quantum creep phenomenon. Extrapolation of the 
theoretical results for dirty superconductors leads to the 
conclusion that the zero temperature magnetic relaxation 
rate scales with the zero temperature normal state con- 
ductivity cr„. Contrary to this, as shown below, the mag- 
netic relaxation rate of Y,i-jPrj^'iBa2Cuj,0'j-s remains 
finite in spite of cr„ (T) ^ at T ^ 0. 

We reported the observation of quantum creep in 
Yi-,j^Pr,j^Ba2CuzO'j^s crystals in a recent Letter [ p3[ . 
In this paper, we present more data and an extended 
analysis of both transport and magnetic relaxation mea- 
surements on two strongly underdoped single crystals of 
Y,4iPr,53Ba2Cu30T-s with Tc « 17 and 21 K, respec- 
tively. The two twinned single crystals were grown by a 



self- flux technique as described elsewhere |^J]. Typical 
dimensions are 0.8 x 0.5 x 0.015 mm^, with the c-axis of 
the crystals oriented along the smallest dimension. 

By performing transport measurements as a function 
of applied current / and magnetic field H, we demon- 
strate the existence of a current independent coupling 
transition temperature T* (H) preceding the crossover to 
quantum creep. The dissipation is ohmic above T* and 
non-ohmic below T* . A picture that arises from these ob- 
servations is that the sample, at T < T*, is divided into 
two macroscopic regions: a layer near the primary face 
(where the current contacts are located) which carries 
most of the transport current, and the rest of the sample 
which remains mostly undisturbed by the current. Inside 
these layers the vortices are coupled, with the correla- 
tion length comparable to the thickness of the respective 
layer. On the other hand, these two macroscopic regions 
are decoupled from each other. As a result, the ratio 
Vp/Vs of the primary to the secondary voltage increases 
with decreasing temperature several orders of magnitude 
over its value in the normal state. 

The crossover to temperature independent (quantum) 
creep takes place at T = Tq{H) < T*{H) and only in the 
top (current-carrying) layer; the rest of the sample con- 
tinues to exhibit thermally activated creep. This leads us 
to the conclusion that the upper layer has a thickness of 
only a few unit cells similar to ultrathin films and multi- 
layers, the only other systems in which quantum creep 
has been observed in transport. Outside this region, 
the vortices have much greater length, which suppresses 
quantum tunneling because the probability of tunneling 
decreases exponentially with the length of the tunneling 
segment. We argue below that the formation of one or 
several very thin channels that carry a current density 
much greater than average may be a result of non-ohmic 
dissipation below T* which can lead to an instability of 
the current distribution with respect to self-channeling. 

In order to confirm that the T-indcpcndcnt dissipa- 
tion is due to quantum tunneling, we performed mag- 
netic relaxation measurements on a similar crystal of 
Yi-xPTxBa2Cu30r-s . The results show a transition 
to T-independent relaxation rate at approximately the 
same temperature Tq[H) as in transport. In addition, 
we were able to determine the characteristic relaxation 
time which characterizes the rate of relaxation uninhib- 
ited by the lack of thermal energy. The value of this 
"escape time" ~ 1 s indicates that the relaxation of the 
magnetic moment is governed by the diffusion of vortices 
from the bulk to the outer edge of the sample. 



II. TRANSPORT MEASUREMENTS 

Transport measurements were performed using the 
"flux transformer" contact configuration (Inset in Fig. ^). 
The current / was injected through the contacts on one 
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FIG. 1. Primary Vp and secondary 14 voltages normalized 
to the total current I and magnetic field H plotted versus 1/T 
for five magnetic fields (0.2, 0.4, 0.6 0.8, and 1 T). The slope 
decreases with increasing field which is parallel to the c-axis. 
Inset: Contact configuration used in our measurements. 

face of the sample and the voltage drops between contacts 
on the same (primary voltage Vp) and the opposite (sec- 
ondary voltage Vs) faces was measured for temperature, 
total current, and magnetic field applied parallel to the 
c-axis of the crystal in the ranges of (1.9 K <T < 20 K), 
(0.3 nA < I < 2 mA), and (0.2 T < H < 9 T) . 
The single crystal, cleaved into a bar-shaped sample, was 
mounted on a single-crystal MgO substrate. {MgO and 
YBa2Cu30'j-s have similar coefficients of thermal ex- 
pansion.) The eight electrodes were fabricated by bond- 
ing 2 mil Au wire to the sample with Ag paste. A typical 
contact resistance was 2 i7 or less. The mean-field su- 
perconducting transition temperature Tco ~ 16.9 K was 
determined from the temperature dependence of the in- 
plane electrical resistance R{T) measured in a magnet- 
ically shielded environment {H < 10^^ G) with a low 
transport current density (< 10 A/crn^) by using the 2D 
Coulomb gas model (2||. At T^q (~ 16.9 K) the resis- 
tance of the sample is 90% that of the normal state value. 



A. Thermally assisted creep 

Figure U gives an overall view of the temperature and 
field dependence of the primary Vp and secondary Vs volt- 
ages. These voltages are normalized to the current and 
field. The convergence of these curves at T x 9 K indi- 
cates a regime where the dissipation is due to the free mo- 
tion of vortices. At lower temperatures, both resistances 
exhibit activated T— dependence with field-dependent ac- 
tivation energies. Due to the resistive anisotropy of the 
crystal, greater current flows near the primary face so 
that Vp > Vs in both the normal and mixed states. At 
even lower T, the primary voltage becomes T indepen- 



dent indicating onset of quantum creep while the sec- 
ondary voltage remains thermally activated. 

At T ^ 9 K, the resistance determined from both the 
primary and secondary voltages is proportional to the 
applied field, i.e., 






(1) 



where Ri = Rp^g = yp,s/I and i?" is the corresponding 
normal state resistance. Equation (|^) describes the free 
flow of vortices near the upper critical fleld Hc2 (T) [^ . 
Fields of 0.2 T and higher are large enough, for this sam- 
ple, to shift the onset of free flux flow regime substantially 
below the zero field Tc « 17 K. 

The system of vortices undergoes a transition into a 
new state at a sharply defined temperature T* {H) which 
is the same for both primary and secondary voltages and 
decreases with increasing H (Figs. 0(a) and p|(b)). For 
T > T* , the activation energies near the primary and 
secondary faces of the crystal {Ep^s = —d\nVp,s/d{l/T)) 
are equal and, therefore, current independent. Below T*, 
both Ep and Es change [both Vp^s{T) curves acquire a 
different slope], but always Es > Ep. The value of T* 
decreases strongly with increasing field. 

It is interesting that the ratio Vp/Vs which changes 
with temperature and field appears to scale with the 
values of T*, as demonstrated in Inset in Fig. ||(a). 
Plotted versus T/T*, the data points for different fields 
form a single curve. Note also that the transition to T- 
independent creep takes place aXT — Tq ^ 0.55 — 0.6 T*. 
Such a strong correlation between Tq and T* is a clear 
indication that the changes in the vortex system which 
occur at T* have a strong impact on the transition to 
quantum creep, or, perhaps, are a prerequisite for such a 
transition. 

Further details of the transformation at T* are pre- 
sented in Fig. 0, where the primary resistance Rp{T) is 
shown for several values of the transport current within a 
range of over two decades (1 /iv4 < I < 250 fiA). The dis- 
sipation at T > T* is ohmic, so that the resistance Rp{T) 
and activation energy Ep are current independent. How- 
ever, the dissipation becomes non-ohmic below T* . At 
low currents, the activation energy below T* is greater 
than above T*, so that the curve Rp{l/T) has downward 
curvature. The activation energy decreases with increas- 
ing current, and, at sufficiently large currents, Ep{I) be- 
comes smaller than it is at T > T* . For large enough 
currents, therefore, Rp{l/T) acquires upward curvature. 
This explains the upward curvature around T* of Rp(T) 
in Fig. E|(a) (large current) and its downward curvature 
in Fig. |[(b) (smaller current). The secondary voltage Vs 
reflects the substantially smaller current density reaching 
the secondary (bottom) surface of the sample. It shows 
some degree of non-ohmicity, but not as pronounced as 
the primary voltage. 



10" 



b I I I I I I I I I I I I I — I I I — n — I I I — n — r 




1000 =; 



0.1 0.2 0.3 0.4 0.5 0.6 

T' (K-1) 
FIG. 2. Arrhenius plots of the primary Vp and secondary 
Vs voltages measured in two different fields and currents; (a) 
0.6 T and / = 100 fiA, and (b) 1.5 T and J = 10 fiA. 
Inset to (a): The "anisotropy" Vp/Vs plotted versus reduced 
temperature T/T* for applied magnetic fields H of 0.2, 0.4, 
0.6, 0.8, 1.0, 1.5, 2.0, 2.5, and 4.0 T. 

The kink at T* in Rp (Figs. |(a), |(b) and |) is similar 
to that observed in four-point resistive measurements on 
Mo77Ge33/Ge and Mo/ Si multilayers j|,^ and oxygen 
deficient YBa2CuzOT-s thin films pO]. However, the 
Mo'jjGe^'i/Ge multilayers exhibit a downward curvature 
in R{T) Pi, while the YBa2Gu307-s films show an up- 
ward curvature pq |; the Mo/ Si shows downward curva- 
ture for some samples and upward curvature for others 
|p|. The data in the Inset in Fig. demonstrate that the 
origin of this contradiction is the current dependence of 
the activation energy below T* . The threshold current at 
which the curvature of R{T) changes sign is material and 
sample specific, which explains seemingly contradictory, 
in this respect, outcomes of different experiments. 

Furthermore, Fig. [^ shows that the crossover temper- 
ature T* is current independent. This, along with the 
fact that T* is the same for both the primary and sec- 
ondary voltages indicates a thermodynamic transition at 
T* rather than a kinetic phenomenon. All these results 
indicate that this thermodynamic transition is between 
a system of decoupled 2D vortices and a system of 3D 
vortices with the macroscopic coherence length along the 
direction of the magnetic field (the c-axis in this case) . 
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FIG. 3. Temperature dependence of the primary resistance 
Vp/I for different values of the total current measured in a 
magnetic field H — 0.2 T. The resistance is ohmic above 
T* and non-ohmic below T*. For clarity, the inset shows the 
same data for two values of current, IfiA and 250/ij4. The 
solid lines are guides to the eye. 

At T > T*, the vortices are not coherent in the di- 
rection of the magnetic field (decoupled) and behave as 
2D "pancakes" . The dissipation mechanism is activated 
hopping of 2D pancakes over potential barriers since the 
activation energies are the same for the primary and sec- 
ondary surfaces in spite of the nonuniform current dis- 
tribution. Figure is a plot of the activation energy 
[/2^ = Ep^s versus field at T > T*. U^°{H) decreases 
monotonically with increasing field. This is expected be- 
cause some of the vortices fill the deepest pinning wells. 
This leads, due to mutual repulsion, to a smoother po- 
tential profile (healing of the random potential) and, on 
average, lower activation energies for the rest of the vor- 
tices, which provide the bulk of the dissipation. The 
decrease of the activation energy due to healing is most 
pronounced at low vortex densities (low H) since the vor- 
tices heal the deepest parts of the pinning profile first. 
The efficiency of this process decreases at higher vortex 
densities (higher H). Correspondingly, the rate of change, 
dU^^ / dH , decreases with increasing field. 

The pancake vortices form coherent lines at T < T* 
and, hence, the activation energy increases {Vp,s curve 
downward). However, the vortex lines do not extend 
through the whole thickness of the sample, so that Vp re- 
mains greater than Vs and the ratio Vp/Vs even increases 
below T* (see Inset in Fig. ||(a)). Moreover, the activa- 
tion energy decreases with increasing current and eventu- 
ally becomes smaller than that for 2D pancake vortices 
(Vp acquires an upward curvature). This is clearly in- 
consistent with the idea that the correlation length Lc is 
limited by "flux cutting" processes |g^. If Lc along the 
c-axis is destroyed by the stress due to the driving 
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FIG. 4. Magnetic field dependence of the activation energy 
[/^^ = dlnVp,s/d{l/T) determined at T > T*. The sohd hne 
is a guide to the eye. 

force, the activation energy would decrease with increas- 
ing current but could not become smaller than it is for 2D 
pancake vortices. It is obvious then, that below T* a new 
channel of relaxation opens up and becomes dominant at 
sufficiently large currents. 

Both the field and current dependences of the activa- 
tion energy at T < T* are consistent with a 3D plas- 
tic creep model based on dislocation mediated motion of 
vortices, similar with diffusion of dislocations in atomic 
solids p§], with the activation energy given by pqj2^: 



Upi{I,B) = U;,{B){l~{^ 



jpl ■ 



(2) 



where B is the magnetic induction, U°i{B) is the acti- 
vation energy at / = and I^'' is the critical current 
corresponding to the plastic motion of the vortices. In 
the limit of small currents, the activation energy U°i for 
the motion of a dislocation in a 3D vortex system can 
be estimated as the energy needed for the formation of a 
double kink over the Peierls barrier i22,E9,pO|; i.e.. 



UZ{B) 



2aoto 



$2 



(%)^ 



^^^iKbB 



(3) 



where Eq = (t-%— )^ is the line tension for a vortex 
aligned along the c-axis, $o = 2.07 x 10^^ Gem? is the 



fluxoid quantum, a^ 



is the vortex lattice con- 



stant, 7 — \ -r^ is the anisotropy parameter, and Xab.c 
is the penetration depth associated with screening cur- 
rents flowing in the ab plane and c-axis, respectively. 

Figure ^a) is a log- log plot of U'^^ versus H, extracted 
from the secondary voltage Vs{T, H) where the current 
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FIG. 5. (a) Field H dependence of the activation energy 
U^^ determined from the secondary voltage as din Vs/dV^/T) 
for T <T*. The solid line is a fit of the data to Eq. @.{h) 
Current I dependence of the activation energy Up^ deter- 
mined from the primary voltage as din Vp/d{l/T) for T < T* 
in a field H = 0.2 T. In the limit of small current, the value of 
Up^ is equal to that of U^^ determined from the secondary 
voltage in the same range of temperature T <T* . The value 
of U is also indicated by the double arrow. 

is very small (/ -^ 0). The data exhibit an H~^/^ de- 
pendence which is characteristic of the motion of a dis- 
location in a 3D vortex structure (Eq. (||)). Figure g(b) 
displays the current dependence of the activation energy 



Up — rj. 



p extracted from the primary voltage Vp{T,I) 



for T <T*, measured in a magnetic field H = 0.2 T. The 
double arrow indicates the value of the current indepen- 
dent t/^^ at T > r* for the same magnetic field. Notice 
that U^^ > U^^ for / < 0.1 mA and U^° < C/^^ for 
/ > 0.1 niA. At this threshold current, Vp{T) changes its 
curvature from downward open to upward open. In sum- 
mary, these results show that the dissipation at T < T* 
is determined by two parallel processes: thermally acti- 
vated motion of correlated vortices (dominant at low cur- 
rents) with the activation energy greater than that for a 
2D vortex, and plastic motion of dislocations (dominant 
at higher currents) with the activation energy smaller 
than that for a 2D vortex. 

A schematic model of the flux flow which transpires 
from these observations is shown in Fig. |6|. The current 
applied through the contacts on the primary face creates 
a nonuniform Lorentz force acting on the vortices. 
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FIG. 6. Schematic representation of the current flow within 
the single crystal. The transport current flows mainly in the 
upper (shaded) layer. The vortices are coupled within each 
layer, but the two layers, current-carrying and "dormant" 
(unshaded) are decoupled from each other. The correlation 
length of vortices in the dormant layer is greater than that in 
the current-carrying layer. 

which is equivalent to the application of a shear stress to a 
fragile solid at T < T*. The shear stress triggers plastic 
flow which is strongest near the primary surface where 
the current density is greatest. On the other hand, the 
vortex lines near the secondary face remains relatively 
undisturbed. This difference of shear stress results in 
asymetric flux growth such that the vortices that grow 
from the secondary surface are longer than those that 
originate on the primary face. These two regions of the 
sample are uncoupled, yet for each of them the vortices 
are coherent over macroscopic distances. In this scenario, 
the large resistive anisotropy measured in the mixed state 
results from the loss of the phase coherence only between 
two macroscopic regions of the sample, not between all 
microscopic layers (such as Cu02 bilayers) as in the 2D 
phase at T > T*. 



B. Quantum creep 

At lower temperatures, another transition at a field 
dependent temperature Tq takes place (see Figs. 1, 2(a), 
and 2(b)). The primary voltage Vp becomes tempera- 
ture independent and scales with the applied magnetic 
field; i.e., the resistance curves Rp = Vp/I normalized 
to the magnetic field Rp/H tend to converge below Tg 
(Fig. IT]). It is important to note that, although the sec- 
ondary voltage (which represents the dissipation in the 
lower section of the sample) does not exhibit the transi- 
tion to T-independent resistance, the activation energy 
is noticeably smaller below Tq (Figs. g(a) and ||(b)). 

Since most of the current flows in a thin layer near the 
primary surface, it is useful to give an estimate of the 
value of the sheet resistance and residual mobility: 
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FIG. 7. Current / dependence of the sheet resistance Ra 
in the quantum creep regime for two values of magnetic field 
0.2 and 0.4 T. 
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Here b « 0.5 mm is the width of the sample and 
£ « 0.3 mm, is the distance between the voltage con- 
tacts. Following the conventional treatment pM of the 
dissipation due to vortex drift, the sheet resistance can 
be expressed in terms of the vortex mobility fires (the 
total mobility of a moving segment, not the mobility per 
unit length): 



f^re 
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(5) 



where c is the speed of light and 4>q is the flux quan- 
tum. For Rn oc H (Fig. |l|), the vortices below Tq 
are characterized by a flnite, fleld independent "resid- 
ual mobility" even in the limit T — > 0. The saturation 
of Vp at low temperatures persists even at the lowest 
current of 0.01 m,A. At low currents, Ra = TZH with 
7^ « 20 mQ/T = 6 X IQ-^h/e"^ T"!; the corresponding 
"residual mobility" fires « 1 x 10^°s/g. 

As shown in Fig. |^, the resistance Ra is current inde- 
pendent at lower currents and increases with increasing 
current at / > 10 fi A. The current dependence of Rn can 
be well fitted with 



Ra=HTZ[l 



1/2 



(6) 



with /o « 20^^. 

These observations clearly indicate that, in strongly 
underdoped Yi^xPfxBa2Cuy,0'j-s , quantum creep be- 
gins to dominate classical, thermally activated creep at 
relatively high temperatures Tq r^ 5 K in H = 0.2 T. A 
factor that may facilitate a transition from classical to 



quantum creep is the high normal state resistivity. Pre- 
viously, we measured the normal state resistivity Pn{T) 
of the same sample by suppressing the superconductivity 
with a large magnetic field ||l^. This showed that the 
normal state of this superconductor is insulating, simi- 
lar to that of PrBa2CuzO'!-s, so that Pn{T) ^ oo as 
T — > 0. The reduced dissipation in the normal core of 
a vortex due to a large normal-state resistivity increases 
the mobility of the vortices ||l^ and, therefore, facilitates 
tunnehng | pO[ |. 

However, current theories cannot be directly applied 
to this system, because they predict the Euclidian action 
Se of the tunneling process to scale to zero with the zero 
temperature normal-state conductivity ||32|; i.e.. 



Se ,, nL, 
n '*e2p„(0)' 



(7) 



where L^ is the length of the tunneling segment. It is 
clear, however, from the small values of the sheet resis- 
tance and residual mobility, Ra oc pres oc exp{— 5b/?1}, 
that the Euclidian action does not tend to zero, but re- 
mains finite at T ^ 0. Hence, the theory of vortex tun- 
neling, as well as the Bardeen-Stephen treatment of vis- 
cosity must be modified for systems which have insulating 
normal state underlying the superconductivity. A serious 
discrepancy between experimental and theoretical values 
of viscosity was also noted in Ref . |1^ by the analysis of 
the relaxation rate in dirty superconductors. 

Previously, the quantum creep in transport measure- 
ments was observed exclusively in thin films whose thick- 
ness did not exceed 30 — 40 A p4|-[l6[. The main rea- 
son for this is the exponential decrease of the probability 
of tunneling with increasing correlation length along the 
field direction. In films, this length is restricted by the 
thickness of the film. Our crystals have a much greater 
thickness, about 1.5 x lO^A. However, the fact that we 
observe a T-independent primary voltage Vp and a ther- 
mally activated Vg indicates that the thickness of the 
"upper" (current-carrying) layer (see Fig. g) is proba- 
bly self-restricted to just a few unit cells, thus facili- 
tating tunneling of such short segments even at rela- 
tively high temperatures (T « 5 K). The correlation 
length is much longer in the rest of the sample (below the 
current-carrying layer) , comparable to the total thickness 
of the sample and, as a result, the vortices do not tunnel 
and the secondary voltage remains thermally activated 
(Figs, m and ||). It is interesting, however, that the tran- 
sition to quantum creep still has an effect on the bulk 
of the crystal, because the activation energy determined 
by the slope of the secondary resistance \dlnRs/d{\/T)\ 
decreases at Tq. 



III. RELAXATION OF MAGNETIZATION 

To verify that the transition to a temperature indepen- 
dent creep, discussed above, is due to quantum tunneling, 
we also performed magnetic relaxation measurements 
on another single crystal of YA7Pr,ziBa2CujOT-s , us- 
ing a SQUID magnetometer over a temperature range 
2 K < T < 2Q K for applied magnetic fields H up 
to 5 T. A small, 3 cm, scanning length was used 
to minimize the variations in field strength inside the 
sample due to spatial inhomogeneities in the magnet 
{5H < 0.048%). The superconducting transition tem- 
perature Tc ~ 21 X of this single crystal was determined 
from the onset of diamagnetism measured in a low mag- 
netic field {H = 10 Oe). The irreversibility temperature 
Tirr for a given H was defined as the temperature above 
which the zero-field-cooled and field-cooled magnetic mo- 
ments are identical. Magnetic relaxation measurements 
were performed by cooling the sample in zero field, ap- 
plying a field H + AH {AH = 0.3 T for all H) parallel to 
the c-axis of the crystal and then reducing it to H. The 
decay of the resultant paramagnetic moment was moni- 
tored for several hours ( w lO'' s ) in constant field H. 
This procedure was used to ensure that the sample was 
in the fully critical state |3^]. The irreversible part of 
the magnetic moment Mirr was obtained approximately 
by subtracting the field-cooled moment from the total 
measured moment. 

From this data we can determine whether the relax- 
ation process also exhibits a transition from thermally 
activated to quantum relaxation at the same temperature 
as the transport resistance. When the relaxation of the 
magnetic moment proceeds as a sequence of uncorrelated 
microscopic steps, each requiring thermal activation over 
an energy barrier, the decay time t^ during which the 
induced moment loses a substantial fraction of its initial 
value can be expressed as: 



Td 



.exp 



(U{H,T) 



T 



(8) 



Here the Boltzman factor reflects the degree of availabil- 
ity of energy U required for an average elementary step 
to proceed and is essentially independent of the physics 
of the relaxation process. The pre-exponential factor 
Tesc is a measure of how rapidly the relaxation would 
proceed, had it not been limited by the unavailability 
of thermal energy. We call Tesc an escape time to dis- 
tinguish it from the microscopic attempt time r^ which 
characterizes the period of vibration of the vortex inside 
a pinning well. The escape time depends on the size of 
the sample and may depend as well on the magnetic field 
and temperature. Factorization of the decay time given 
by Eq. (0) is meaningful as long as the Boltzman factor 
exp(C//T) ^ 1, so that it dominates the temperature and 
field dependence of Td- Since the activation energy van- 



ishes near Tc, the definition of Td can be specified further 
by taking a hnear T— dependence of the effective barrier: 



U{H,T) 






(9) 



where T^^ is the temperature at which the effective acti- 
vation energy vanishes. It is commonly taken to be equal 
to the critical temperature, but may be smaller than T^ 
and close to the irreversibility temperature. Thus, the 
decay time has the form: 



Td = Tesc exp <{ [/o(-ff) ( - - — 



(10) 



We want to emphasize that Eq. (ph is more general 
than any particular dynamic model of the relaxation pro- 
cess driven by fiuctuations. Therefore, it can also be ob- 
tained within the commonly used model in which the re- 
laxation is described as a decay of the average supercur- 
rent J determined by current-dependent activation en- 
ergyC/(J,B,T): 



It 



-K exp 



U{J,H,T) 
T 



(11) 



Since U{J) increases with decreasing current, this equa- 
tion can be integrated by the method of steepest descent: 



Kt 



^J'exp/^ 
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\dU/dJ\ 
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(12) 



All unknown parameters can be absorbed into one t^sc 
so that an approximate solution of Eq. (O) has the form 



T 



= ln 



t 



(13) 



from which Eq. (bI) immediately follows. Equation (13) 

was obtained earlier |34|] by a more circuitous derivation. 

With U{J) given by the collective creep model [ p2| , 



U{J,H,T)^ 



U{H,T) 



^ -1 



Eq. (O) gives the following time dependence of J: 



J{t)^jJl + i^^\n{t/Tes 



-\lv 



(14) 



(15) 



On the other hand, for an arbitrary C/(J) in Eq. (|13|), the 
initial decay (^Jc — J <IC Jc) is linear in logarithm of time 
(Kim- Anderson formula): 



^(O-^c(l-(^)ln — 



(16) 
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FIG. 8. Time t dependence of the relaxation of the irre- 
versible magnetic moment Mirr measured at different tem- 
peratures in a magnetic field H — 0.8 T. To avoid clutter, 
only a few representative curves are shown. 

where U = Jc\dU/dJc\, the critical current is determined 
by the condition U{Jc) — 0, and it is assumed that 
dU /dJ ~ const. dX J — Jc- 

Due to the slowness of relaxation, the decay time Td 
cannot be directly determined by monitoring the relax- 
ation of the induced moment until it loses a substantial 
fraction of the initial value. An alternative method is 
to estimate the decay time by extrapolating the initial 
decay of J(i) to lower current values. Specifically, when 
the initial decay is described by Eq. (|l6|), we define Td 
from the condition J{Td) = 0. Comparing this definition 
with the collective creep formula, Eq. (15), we see that 
the so defined Td corresponds to a decay to the level of 
Jc/{^ + J^)^'". The current density J can be experimen- 
tally determined by the irreversible part of the magnetic 
moment Mirr{t) oc J{t). 

Representative semilog plots of AIirr{t) oc J{t) as a 
function of time t for several temperatures measured in a 
field H ~ 0.8 T are shown in Fig. |[ Within a decade of 
time lO'' — 10'* s, the relaxation curves can be well fitted 
to: 



Mirr = a — 61n(t/to), 



(17) 



where to is an arbitrary unit of time. The decay time for 
which Mirr{Td) = is then given by: 



Td = to exp I - 1 



(18) 



With this definition, Td is universal and does not depend 
on the choice of to. 

Figure ^ shows the decay time calculated according to 
Eq. (nS) and plotted against the inverse temperature for 
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FIG. 9. Decay time Td versus 1/T for several values of mag- 
netic field {H = 0.1, 0.2, 0.6, 0.8,and 1.2 T). The slope de- 
creases with increasing field. The straight line extrapolations 
of the Arrhenius type dependence converge at Tcr ~ 19 K and 
Td = Tesc ~ 1 s. The saturation of Td at the level 10^^ — lO^^s 
is due to quantum creep. 

different values of magnetic field. A comparison of these 
data with Td{T) given by Eq. (10) allows us to deter- 
mine whether the relaxation crosses over from activated 
to non-activated dependence, and at what temperature. 
At higher temperatures, indeed, the data display an Ar- 
rhenius dependence with a slope d\nTd/d{l/T) decreas- 
ing with increasing field. This trend is consistent with 
the field dependence of the activation energy in trans- 
port measurements. Fig. |[ 

It is important to note that the values of the activation 
energy determined by the slopes d\iiTd/d{l/T) are not 
universal and depend on the criterion used to define the 
decay time (see Eq. (p^). However, the pre-exponential 
factor Tesc is universal and can be determined by the 
extrapolation of the Arrhenius dependence of Td{T) to 
the temperature T^r- Indeed, a linear extrapolation of 
the data in the activated region to higher temperatures 
(straight lines in Fig. 0) shows that the lines converge 
at Tcr ~ 19 K (which is consistent with the value of 
Tc or Tirr)- The point of convergence corresponds to 
Td = Tcsc ~ 1 s. This is an extremely large character- 
istic time of relaxation in comparison with the attempt 
time which is typically assumed to be of the order of 
10~^ — 10""'^^ s. The value of the escape time can be esti- 
mated from the following consideration. The decay of the 
induced moment occurs when vortices leave the sample 
|p5| . Over long times (in comparison with the micro- 
scopic time scale) any type of rearrangement of vortices 
reduces to diffusion. Correspondingly, the escape time 
can be estimated as the time required for a vortex to 
diffuse from the bulk to the outer edge of the sample: 









Wm* 



(19) 



where R is the characteristic size of the sample in the 
direction of diffusion (in the a — b plane in our case), 
and r>i, is the diffusion coefffcient determined by the 
attempt frequency uja and the average elementary vor- 
tex hopping distance ia- With Tcsc ~ 1 s and R^ ~ 
10~^ — lO"'^ cm'^ (for the crystal we measured), Eq. ( |l9| ) 
gives Vy - 10^2 - 10-3 cm'^/s. This value of V^ is 
consistent with an elementary step of the order of the 
correlation length 4 ~ 100 A and uja -- lO^" - 10^ 5"^ 
We define the effective mass m* of a segment of the vor- 
tex line through the uncertainty principle, tOa ~ h/m*£^. 
With these estimates, the effective mass of the diffusing 
vortex segment is 10^ — 10'^ times of the electron mass. 



A. Quantum creep 

At lower temperatures, the decay time saturates 
at a roughly temperature and field-independent level 
(Fig. 1^). The crossover temperatures Tq{H) from trans- 
port (Fig. |l|) and magnetic relaxation (Fig. ||) measure- 
ments are very close in spite of a very large difference in 
the currents involved in these measurements. The fact 
that the transition to a temperature independent dissi- 
pation takes place in both transport and magnetic relax- 
ation processes, and at approximately the same tempera- 
ture in a given field indicates that both phenomena have 
a common origin. 

In the regime of quantum relaxation, the relaxation 
rate is limited by the probability of tunneling as deter- 
mined by the Euclidian action Se- Similar to Eq. (pi]), 
the relaxation rate of the supercurrent can be expressed 
in terms of the current dependent SsiJ)' 



dJ J. ( Se{J,H) 



which has a solution similar to Eq. ([l^): 



Se{J) 

n 



In 



(20) 



(21) 



Provided that dSs/dJ — const at J^, where Jc is deter- 
mined by the condition Se{Jc) = 0, the initial decay has 
same linear in logarithm of time dependence as in the 
classical case (see Eq. (16)); i.e.. 



Jit) ^Jcil 



h , t 

— In 

OO Tcsc 



(22) 



where So = Jc\dSE/dJc\- The decay time determined by 
the extrapolation of Eq. (03) to J{Td) = is given by: 



Td 



• CSC exp <; — 



^0 



(23) 



The value of the escape time should be similar to that 
in the classical regime since it is determined by the dif- 
fusion uninhibited neither by the lack of the thermal en- 
ergy, nor by the small tunneling probability. Therefore, 
we can estimate 5*0 from the data of Fig. S: 






Td 



25. 



(24) 



This value is comparable, but somewhat smaller than 
those reported for other systems |l^ . 



IV. SUMMARY AND SPECULATIONS 

We have observed quantum creep in underdoped 
Yi-xPrxBa2Cuj,0'j-s crystals using both transport and 
magnetic relaxation measurements. The transition to 
quantum creep is preceded by a coupling transition which 
leads to non-ohmic dissipation. The evidence presented 
in previous sections lead us to a picture of the current 
density distribution shown in Fig. g. Most of the trans- 
port current is confined to a very thin layer below the 
current contacts. This current-carrying layer is decou- 
pled from the rest of the crystal, where the vortices are 
mostly undisturbed by the current and are coherent over 
a macroscopically long distance, perhaps comparable to 
the thickness of this "dormant" layer which is practically 
the same as the thickness of the sample. 

The transition from thermally activated to temper- 
ature independent dissipation takes place only in the 
current-carrying layer. The relatively large tunneling 
probability which makes possible the observation of this 
crossover at T w 5 K in these crystals is due to very short 
tunneling segments, large normal state resistivity, and a 
large current density. In the rest of the sample, the vor- 
tices are much longer and their tunneling is suppressed as 
manifested by activated T-dependence of the secondary 
voltage down to the lowest temperature. This is con- 
sistent with the fact that quantum creep was previously 
observed only in transport experiments on ultrathin films 
and multilayers with a thickness no more than 30 — 40A. 

Magnetic relaxation measurements substantiate that 
this T-independent resistance is due to quantum creep. 
The decay time of the magnetic moment becomes T- 
independent at approximately the same temperatures (in 
a given magnetic field) as in transport (see Figs. |l|, |[ and 
^ . We also determined the characteristic relaxation time 
T^sc which turns out to be very large ~ 1 s in comparison 
with the microscopic attempt time. 

While the finding of nonvanishing resistance at T ^ 
in a crystal and its correlation with nonvanishing mag- 
netic relaxation rate is important and has never been 
observed before, our results raise also another impor- 
tant question. In the normal and mixed state above 
the temperature of the coupling transition T* , the crys- 
tals of Yi^xPrxBa2Cu^0-j-^s are not very anisotropic. 



In the sample with the length L k, \ mm and thick- 
ness D « 0.015 m?7i, the ratio Vp/Vs ~ 2, (see Inset in 
Fig. 11(a)) so that the transport current fills fairly uni- 
formly the whole cross-section. On the other hand, the 
transition to quantum creep indicates that below Tq the 
current-carrying volume collapses into a thin layer, pos- 
sibly just a few unit cells thick. This favors the quantum 
creep for two reasons: minimum length of the vortex seg- 
ments and maximum current density which reduces the 
height of pinning barriers. This opens the question of the 
nature of such a drastic self-channeling of the transport 
current. 

A possible answer to this question is related to the non- 
ohmic, current-dependent resistive anisotropy. This can 
be illustrated by the following qualitative dimensional 
considerations. According to local electrodynamics, in 
the sample with thickness D (Inset in Fig. P, the trans- 
port current mostly flows within a layer of thickness £'e// 



D. 
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— ; V 

V 



L \PabJ 



(25) 



If the effective anisotropy increases with increasing cur- 
rent density j then the transport current exhibits a ten- 
dency to channel itself into an increasingly narrow layer 
as follows. Let us assume that to the lowest order in j, 
the anisotropy can be written as: 



JO, 



where 



] 



D, 



if 



(26) 



(27) 



and / is the total current. Combining Eqs. ( p5| ) - (p7[), we 
obtain the following equation for the effective thickness 
of the current distribution: 



Deff^ 



D 



r?o(l + ip^)' 



(28) 



This equation has the solution 



Deff = 



D 

2^ 




P 



(29) 



where Ims — Djo/'^Vo- ^^ the total current / increases, 
Deff (I) gradually decreases until it reaches half of its zero 
current value D/rjQ at / = I^ns- For / > lins, Eq. (pS) 
does not have a real solution except for I?e// = . Thus, 
lins is a threshold of instability (the corresponding value 
of Tjins = riilins) = 277o). For / > I,ns there is no sta- 
ble current distribution with macroscopic thickness. The 
current-carrying layer compresses itself until it is a few 



10 



unit cells thick, or until the current density approaches 
the critical value. The existense of such an instability 
would have significant implications to our understanding 
of the electrical transport in layered superconductors and 
their applications. 



V. ACKNOWLEDGMENTS 

This research was supported at KSU by the National 
Science Foundation under Grant Nos. DMR-9601839 and 
DMR-9801990, and at UCSD by U.S. Department of En- 
ergy under Grant No. DE-FG03-86ER-45230. 



[1 

[2; 

[3] 
[4] 

[6; 



[7: 

[9 

[lo; 



[11 

[12; 
[is; 

[14 

[is; 

[16 

[ir 
[is; 



I. Giaever, Phys. Rev. Lett. 15, 825 (1965). 

W. R. White, A Kapitulnik, and M. R. Beasley, Phys. 

Rev. Lett. 66, 2826 (1991). 

N. Y. Fogel, V. G. Cherkasova, O. A. Koretzkaya, and A. 

S. Sidorenko, Phys. Rev. B 55, 85 (1997). 

H. Safar, P. L. Gammel, D. A. Huse, S. N. Majumdar, L. 

F. Schneemeyer, D. J. Bishop, D. Lopez, G. Nieva, and 

F. de la Cruz, Phys. Rev. Lett. 72, 1272 (1994). 

R. Busch, G. Ries, H. Werthner, G. Kreiselmeyer, and 

G. Saemann-Ischecnko, Phys. Rev. Lett. 69, 522 (1992). 
H. Safar, E. Rodriguez, F. de la Cruz, P. L. Gammel, 
L. F. Schneemeyer, and D. J. Bishop, Phys. Rev. B 46, 
14238 (1992). 

R. A. Doyle et al. Phys. Rev. Lett. 77, 1155 (1996). 

C. D. Keener et al. Phys. Rev. B 55, R708 (1997). 

Y. Eltsev, W. Holm, and O. Rapp, Phys. Rev. B 49, 
12333 (1994). 

A. C. Mota et al. Progress in High Temperature Super- 
conductivity, Proceedings of ICTPS 90, ed. by R. Nicol- 
sky (World Scientific, Singapore, 1990), VI. 25. 
A. C. Mota et al. Physica C, 185-189, 343 (1991) and 
references therein. 

A. J. J. van Dalen, R. Griessen, S. Libbrecht, Y. 
Bruynseraede, and E. Osquiguil, Phys. Rev. B 54, 1366 
(1996); and references therein. 

A. F. Th. Hoekstra, R. Griessen. A. M. Testa, J. cl Fat- 
tahi, M. Brinkmann, K. Westerholt, W. K. Kwok, and 
G. W. Crabtree, Phys. Rev. Lett. 80, 4293 (1998). 
Y. Liu, D. B. Haviland, L. I. Glazman, and A. M. Gold- 
man, Phys. Rev. Lett. 68, 2224 (1992). 

D. Ephron, A. Yazdani, A. Kapitulnik, and M. R. 
Beasley, Phys. Rev. Lett. 76, 1529 (1996). 

J. A. Chervenak and J. M. Valles, Jr., Phys. Rev. B, 54, 

R15649 (1996). 

A. Gerber and J. M. Franse, Phys. Rev. Lett. 71, 1895 

(1993). 

G. A. Levin, T. Stein, C. C. Almasan, S. H. Han, D. A. 

Gajewski, and M. B. Maple, Phys. Rev. Lett. 80, 841 

(1998). 



[19 

[2o; 

[21 

[22; 

[23; 
[24; 
[25; 

[26 
[27 
[28; 
[29; 

[30 
[31 
[32; 
[33 
[34; 
[35; 
[36 



J. Bardeen and M. J. Stephen, Phys. Rev. 140, 1197A 
(1965). 

A. O. Caldeira and A. J. Leggett, Ann. Phys. (N.Y.) 149, 
374 (1983). 

G. Blatter, V. B. Geshkenbein, and V. M. Vinokur, Phys. 
Rev. Lett. 66, 3297 (1991). 

G. Blatter, M. V. FeigeI'man, V. B. Geshkenbein, A. L 
Larkin, and V. M. Vinokur, Rev. Mod. Phys. 66, 1125 
(1994). 

T. Stein, G. A. Levin, C. C. Almasan, D. A. Gajewski, 
and M. B. Maple, Phys. Rev. Lett. 82, 2955 (1999). 
L. M. Paulius, B. W. Lee, M. B. Maple, and P. K. Tsai, 
Physica C 230, 255 (1994). 

P. Minnhagen and P. Olsson, Phys. Scripta T 42, 9 
(1992). 

X. G. Qiu, B. Wuyts, M. Maenhoudt, V. V. Moshchalkov, 
and Y. Bruynseraede, Phys. Rev. B 52, 559 (1995). 
D. Lopez, G. Nieva, F. de la Cruz, H. J. Jensen, and D. 
O'Kane, Phys. Rev. B 50, 9684 (1994). 
J. P. Hirth and J. Lothe, Theory of Dislocations (John 
Wiley & Sons, New York, 1982), Chap. 15. 
Y. Abulafia, A. Shaulov, Y. Wolfus, R. Prozorov, L. 
Burlachkov, Y. Yeshurun, D. Majer, E. Zeldov, H. Wuhl, 
V. B. Geshkenbein, and V. M. Vinokur, Phys. Rev. Lett. 
77, 1596 (1996). 

V. B. Geshkenbein, A. L Larkin, M. V. FeigeI'man, and 
M. V. Vinokur, Physica C 162-164, 239 (1989). 
M. Tinkham, Introduction to Superconductivity 
(McGraw-Hill, Inc., 1996), p. 166. 

G. Blatter and V. B. Geshkenbein, Phys. Rev. B 47, 2725 
(1993). 

For a review see Y. Yeshurun, A. P. Malozcmofi', and A. 
Shaulov, Rev. Mod. Phys. 68, 911 (1996). 
V. B. Geshkenbein and A. I. Larkin, Zh. Eksp. Teor. 
Phys. 95,1108 (1989); Sov. Phys. JETP 68, 639 (1989). 
M. V. FeigeI'man, V. B. Geshkenbein, and V. M. Vi- 
nokur, Phys. Rev. B 43, 6263 (1991). 
G. A. Levin, J. Appl. Phys. 81, 714 (1997) 



11 



